Nonuniversal Quasi-Long Range Order in the Glassy Phase of Impure 

Superconductors 



On 
On 
On 



3 

(N 



43 

o 



CO 



I 

O 

o 



> 
in 



o 

On 
c3 



-a 
c 

o 
o 



X 
S3 



Thorsten Emig, Simon Bogncr and Thomas Nattermann 
Institut fur Theoretische Physik, Universitat zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 

(February 1, 2008) 

The structural correlation functions of a weakly disordered Abrikosov lattice are calculated for 
the first time in a systematic RG-expansion. In the asymptotic limit the Abrikosov lattice exhibits 
still quasi long range translational order described by a nonuniversal exponent rj G which depends 
on the ratio of the renormalized elastic constants re = cee/E-u of the flux line (FL) lattice. Our 
calculations show clearly three distinct scaling regimes corresponding to the Larkin, the manifold and 
the asymptotic Bragg glass regime. The manifold roughness exponent £ rm (re) is also nonuniversal. 
Our results are at variance with those of the variational treatment with replica symmetry breaking 
which allows in principle an experimental discrimination between the two approaches. 



PACS numbers: 74.60. Ge, 05.20.-y 

Conventional type-II superconductors show in addition 
to the flux repulsing Meissner state a second supercon- 
ducting (Abrikosov) phase in which the magnetic induc- 
tion B enters the material in the form of quantized flux 
lines (FLs) which form a triangular lattice. Each FL car- 
ries a unit flux quantum (j>Q = hc/2e. The Abrikosov 
lattice is characterized by a non-zero shear modulus cqq, 
which vanishes at the upper and lower critical fields, H C2 
and H Cl , where continuous transitions to the normal and 
the Meissner state, respectively, occur. Abrikosov in his 
mean-field solution treats FLs as stiff rods. 

Thermal fluctuations roughen the FLs resulting in a 
melting of the Abrikosov lattice close to H C1 and H C2 , re- 
spectively, because of the softening of cqq . This applies in 
particular to high-T c materials with their elevated transi- 
tion temperatures and their pronounced layer structures 
0. At present, it is not clear, whether the transition to 
the normal phase at high field happens in these materials 
via one or two transitions. However, melting of the FL 
lattice has clearly been observed experimentally . 

It is well known, that in addition to thermal fluc- 
tuations in type-II superconductors the effect of disor- 
der has to be taken into account since FLs have to be 
pinned in order to prevent dissipation from their motion 
under the influence of an external current. Randomly 
distributed pinning centers lead indeed to a destruction 
of the Abrikosov lattice, as has been first demonstrated 
by Larkin using perturbation theory Q. In particu- 
lar, he found an exponential decay of the correlations 
of the order parameter for translational long range order 
^g(x) = e lGu ( x ) on length scales larger than a disor- 
der dependent Larkin length L^. Here G and u denote 
a reciprocal lattice vectors and the displacement field 
of the FL lattice, respectively, and x = (xj_,z) is a d- 
dimensional position vector. 

However, as was first shown by Nattermann ||, in 
treating the interaction between the FL lattice and the 
disorder, it is crucial to keep the periodicity of this inter- 



action under the transformation u — > u + R, where R is a 
vector of the Abrikosov lattice. This symmetry, which is 
abandoned in perturbation theory [Q and in the so-called 
manifold models ||, leads to a much slower, logarithmic 
increase of the elastic distortions with the system size 
L §J||. This results in a power law decay of the pair 
correlation function Cg(x) = ($gW*-g(0))' In par- 
ticular, Giamarchi and Le Doussal || calculated Cg(x) 
using (i) a variational treatment for the triangular FL 
lattice and (ii) a functional renormalization group (FRG) 
for a simplified model using a scalar displacement field 
u only. In both cases they found Cg(x_l,0) ~ |x^|~'' G 
with fjG — A(A — d) where d denotes the space dimen- 
sionality. Go is one of the smallest reciprocal lattice 
with G a = 4vr/V3, a = (2<j) /\/3B) 1/2 denotes the lat- 
tice spacing and A is a universal constant (A = 1 and 
A = 7r 2 /9 « 1.1 for the treatment (i) and (ii), respec- 
tively). Thus, even with (weak) disorder there is a quasi- 
long range ordered flux phase which shows Bragg peaks, 
the " Bragg-glass" ||. This result is believed to be valid 
in an extended region of the H — T phase diagram up to 
a phase boundary, where the occurrence of unbounded 
dislocation loops - ignored so far - leads to an instability 
of the Bragg glass phase |P-[h4| . 

The power law decay of Cg(x) is reminiscent of the 
situation in pure 2D-crystals where in the solid phase 
Cg(xj_) ~ |x^|~' ?G . This solid phase is indeed a line 
of critical points with ijg = TG 2 (1 + re _1 )/(47rcii) and 
K = cqs/ch- The cu represent the renormalized elas- 
tic constants which have a finite temperature dependent 
value. At the melting temperature T m rjo, = g(l ~ ^ 2 )- 
p"5| . Similarly, the 2D disordered crystal close to the 
glass temperature T g (i.e. r/G < 2) is described by a line 
of fixed points which depends on re ph[ | . 

On the contrary, the Bragg glass phase is characterized 
by a zero temperature fixed point where the ratio between 
the disorder strength and the clastic energy vanishes on 
large scales L as A*L rf ~ 4 . It has been suggested in M, 
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that the fixed point value A* might be universal resulting 
in a universal coefficient A of rj G . In the present treat- 
ment we will show however, that a FRG treatment of 
the full triangular model leads indeed to a nonuniversal 
K-dependent value of rj G — A*(/i)(Ga) 2 . The situation 
is therefore qualitatively similar to that of 2D crystals 
at the melting temperature. Quantitatively, the effect 
reflects the contributions from the interaction between 
different Fourier modes, which are not considered in the 
variational treatment ]7]||. With n dependent in gen- 
eral on B and T, the observation of a ficld-dcpcndcnt 
rjG would yield the opportunity to judge the validity of 
different approximation schemes under debate p7| . 

Since in the Bragg-glass phase dislocations in the vor- 
tex lattice can be neglected |]9|-|i"4| , its configurations are 
described completely by the elastic displacement field 
u(x). The vortex density is then given by p(x, u) = 
J2u 3{ x ± — R — u(R, z)). The impurities are modeled by 
a Gaussian random potential U (x) with two-point corre- 
lation C/(x)C/(0) = ^ 2 9^±/0S(z). Here 7 = / p 2 in n^ 2 , 
/pin denotes the individual pinning force, rn the impurity 
density, £ the maximum of the coherence and disorder 
correlation length JlJ and g is a short ranged function. 
On scales larger than the London penetration depth A 
the FL lattice can be described by local elasticity theory, 
leading to the Hamiltonian 

U = ~ J d d - 2 zdW {c u (Vi • u) 2 + c 66 (V ± x u) 2 + 

+ c 44 (V z u) 2 }+ J d d xE pin (u,x), (1) 

where x^ = (x, y), and the z component has been gener- 
alized to a d — 2-dimensional coordinate to allow for an 
e-expansion about 4 dimensions. The random potential is 
defined as E p i n (u,x) = J7(x)p(x, u). The Hamiltonian is 
invariant under simultaneous mapping of shear onto com- 
pressional modes (rotation of u by 90°) and permutation 
of the corresponding moduli (k — > For E p - ln (u, x) 

this invariance holds in a statistical sense. Thus the parti- 
tion function and all thermodynamic quantities also show 
this symmetry. We therefore can restrict the calculations 
to < k < 1 in the following. Results for k > 1 corre- 
spond to those for k^ 1 with u rotated. 

Performing the disorder average and discarding 
spatially rapidly oscillating terms ||, one obtains 
for the disorder correlator, which is defined by 
E p - ln (u, x)iJpi n (0, 0) = i?(u)(5(x), a periodic function 
with the Fourier representation 



,B 



i?(u)=7e 2 - T ^.g(Ge)cos(Gu) 



(2) 



where B/(j)o is the mean area density of vortices and 
g(G£) is the Fourier transform of g(x±/£). 

To obtain the disorder averaged configuration of the 
FL lattice on a particular length scale, one has to take 
into account the renormalization of i?(u) by fluctuations 
on shorter length scales. This can be done systemati- 
cally by a functional renormalization group (FRG) for 



the replica Hamiltonian resulting from Eq. ([!]) after the 
disorder average. Because of the statistical invariance of 
the the replica Hamiltonian with the inter-replica cou- 
pling (|^) under a shift of u by an arbitrary vector field 
inducing a compression, shear and/or tilt of the FL lat- 
tice, there is no renormalization of the elastic moduli Jl8| . 
Therefore the temperature obeys the exact flow equation 
dT/dl = (d - 2)T leading to a T = fixed point for 
d > 2. Notice however that in the original model of (|l|) 
the statistical invariance is not exactly fulfilled on length 
scales smaller than a, leading to a small renormalization 
cu — > cu of the elastic constants which will be consid- 
ered from now on as effective parameters. In our FRG 
only coordinates are rescaled as x — > exp(d£)x to keep 
the cutoff A fixed with dl the infinitesimal width of the 
momentum shell. Because of the dispersion of the elastic 
constants on scales smaller than the penetration depth 
A, we have to choose here A ss 2ir/\. Fluctuations on 
smaller scales can be ignored if the Larkin length (see 
below) is much larger then A, i.e. for weak disorder. For 
larger disorder one has to take into account the disper- 
sion cu and C44 which will result in a more complicated 
cross-over but not affect the asymptotic behavior. The 
flow equation for R(u) can then be derived along the 
lines of Refs. ]l9| , p0] . Contrary to previous cases, after 
suitable rescaling of R(u), we obtain 



dR(u) 



= eR(u) 



a" 
~2 



(d 2 x R) 2 + {d 2 y R) 2 + 2{d x d y R) 2 



dl 

2A(d 2 x + d 2 y )R- u - (d 2 x R-d 2 y RY+A{d x d y Rf \ (3) 



with the dimensionless parameter A = —d x d x R(0) = 
~d y d y R(0) and d x = d/du x etc.. The last equality as 
well as d x d y R(0) = follow from the requirement of 
hexagonal symmetry for R(u). Performing the momen- 
tum shell integrals to lowest order in e, one obtains the 
anisotropy parameter 
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21n( K ) 



(4) 



i.e., < S < 1 for any ratio n — cqq/c\\. In the isotropic 
case Cu = C66 (S = 0) the flow Eq. (||) reduces to that of 
Ref. po|| , if cu — > 00, as often assumed for FL lattices, 
we have to put (5=1. 

From Eq. (|) we find that the bare, unrenormalized 
dimensionless parameter Aq is given by 



A « 10- 3 A d - 4 7 ^ 



1 



C44C66 



-^(l + G(e/a)), (5) 



where we have assumed g(G£) = £ 2 0(1 - G£) with 9 the 
step function and £ <C a to evaluate the sum over G. The 
Larkin length (parallel to B) in three dimensions is given 
by w ^ 2 /A Aa 2 jjj. The condition of weak disorder 
mentioned above reads then A a 2 <§C £ 2 . 

Now we integrate Eq. (||) to obtain the renormalized 
function i?(u) on all length scales L = exp(Z)/A including 
the fixed point R*(u) for L — > 00. With the bare corre- 
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FIG. 1. RG flow of A(/) through three regimes obtained by 
numerical integration of Eq. 

lator of Eq. (||) showing the full symmetry of the triangu- 
lar lattice - translation, 6-fold rotational axis, 3 mirror 
lines - and the flow of Eq. (||) preserving these sym- 
metries as it ought to, the set of possible solutions is 
restricted to functions with the full lattice symmetry on 
every length scale. Even for i?*(u) an analytic solution 
not being very obvious, our approach is numerical and 
straightforward: The Ansatz assuring translational plus 
inversion invariance is R(u) = J2g Bg cos(Gu) with the 
sum running over the reciprocal lattice. An infinite set 
of coupled nonlinear but ordinary flow equations for the 
coefficients Rg is derived. Rather than solving directly 
for the fixed point i?Q , the flow equations are integrated 
numerically with the remaining point group symmetries 
exploited here. Convergence to a fixed point from a 
huge bassin of attraction is observed. Of special inter- 
est is the flow of A(7) since it determines the scale de- 
pendence of the renormalized propagator A c s(q)/q 4: with 
A e g(q) = A(ln(A/q))(q/A) 4 ~ d and thus the exponents 
fjc and Crm- Obviously, the numerical integration has 
to be restricted to a finite set of Rg with |G| < |G C |, 
but high accuracy of the flow of A(l) can be restored 
by including up to 360 different coefficients and a fi- 
nite size scaling like extrapolation. The accuracy for 
L — > oo is controlled via the exact (^-independent) re- 
lation ei?(0) = (aA) 2 stemming from (||). 

Two different but typical kinds of behavior of A (I) are 
shown in Fig. [j]. Depending on the ratio £/a, two or 
three scaling regimes with different roughness exponent 



£ defined by A(Z) ~ exp(2(7) can be clearly identified. 
For £/a -C 1, one finds a first crossover from the ran- 
dom force (RF) regime (( = e/2) to the random man- 
ifold (RM) regime (C = Crm) at the Larkin length 
where u ~ £ and a second one at the positional correla- 
tion length L a w L^a/ S y ) 1 ^ rm where u ~ a to the Bragg 
glass regime (( — 0). If £/a ~ 1, L a — > and the 
RM regime disappears. Just like in the 1-d (or CDW) 
case |y a cusp in the second derivatives of R(u) devel- 
ops leading to the fixed point shape shown in Fig. |^. 
The length scale L c on which the cusp appears can be 
determined analytically from the divergence in the RG 
flow of d%R(0) = <9 4 i?(0) = 3d£d%R(0). As can be ex- 
pected, L c is related to the Larkin length known from 
perturbation theory by a numerical factor only. But in- 
terestingly we obtain the typical ratio L c /L^ w 100, i.e., 
the cusp appears clearly beyond the Larkin scale. The 
nonwaiver sality of the exponents Crm(' t ) an d Vg( k ) can 
be determined from A(l) via the slope in the RM regime 
and A*, respectively. Their dependence on k is presented 
in Fig. | 




FIG. 2. Characteristic of the Bragg glass: a cusp in 
dx ± R*(u) at the lattice sites. 

With the numerical value for A*(rt) at hands 
we can now calculate the displacement correlations 
B a b(x) = ([«o(x) - u a (0)][u b (x) - u b (0)]) in the Bragg 
glass phase. Introducing the rescaled z-coordinates z t = 
(c*66/C44) ly/2 Z, zi = {cwjc^Y^z and h(t) — t~ 2 ln(l + 
t 2 ), we find for x 3> L a 
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and B yy (x) follows from B xx (x) by permuting en and 

C66- These correlations lead to the translational order with fjo, = A*(k)(oG) 2 and the geometrical factor 
correlation function Cg(x) = (exp(£G[u(x) — u(0)])), 

which reads r A */ w s~<\2 / i- 

A*( K )(aG) 2 ( ^ 2 1 



C G (x) ~ g G LT(x\ + z 2 )-^/^ 1+ ^{x\ + z 2 ) -^/(2(i+i/k)) 

(8) 
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which describes completely the angular dependencies of 
the translational order. Note that the factor gc goes 
to 1 in the limit z — > oo. Therefore, in this limit the 
dependence of Cg(x) on the reciprocal lattice vector G 
remains only in the exponent f/G. Moreover, it is inter- 
esting to note that the exponents of the algebraic decay 
in Eq. (||) depend on the elastic moduli as soon as z is 
finite even without taking into account the nonuniversal- 
ity of the exponent fjc itself. Neglecting the non-trivial 
renormalization of A* (re), in the case z = the above 
formulae reduce to those found in ftl. 




FIG. 3. Variation of the exponents £ rm and r/a with re. 

The re-dependence of ^g and £ rm (Fig 3), the 
correlation functions B a b(x) (Eqs.(6),(7)) and Cg(x) 
(Eqs.(8),(9)) and the cross-over from the Larkin via the 
random manifold to the Bragg glass regime (Fig. 1) are 
the main results of this paper. 

In isotropic superconductors at low temperatures, 
where flux lines interact via central forces, on has < 
k < 1/3. re ~ 1/3 for A < a, i.e. for fields close to H Cl , 
and k -t for H —* H C2 . For most of the field region 
re (j>o/16Tr\ 2 B,. Thus, an increase of the external field 
from H Cl to H c . 2 should result in an increase of ?7g and 
a decrease of £ rm . Numerically, the effect is small, since 
r^Go ranges from 1.143 to 1.159 and ( rm from 0.1737 to 
0.1763 in this re range. Thus it will be probably hard 
to detect this effect. At higher temperatures, where the 
flux line interaction is considerable influenced by thermal 
interaction, as well as in anisotropic superconductors, 
the above inequality for re may not be longer fulfilled. 
Clearly, in the latter case also our starting Hamiltonian 
(1) would have to be modified. 

Contact to the neutron scattering experiment is made 
by S(G + q) = Jd 3 cce iqx CG(x), the structure factor, 
which is proportional to the scattered intensity. For the 
marginal cases re = and re = 1 we find from Eq. (0) 



i.e., <S(G+q) exhibits Bragg peaks. For re > the Fourier 
transform can be easily done numerically. 

To summarize, we have shown that contrary to pre- 
vious claims, the translational quasi-long range order in 
the Bragg glass phase of impure type-II superconductors 
is described by a nonuniversal power-like decay of the 
order parameter correlations. In particular, the decay- 
exponent fjc depends on the ratio re = cqq/cu of the elas- 
tic constants, similar to 2D pure crystals at their melt- 
ing temperature. For weak disorder we find a crossover of 
the structural correlation functions from a Larkin-regime, 
where perturbation theory applies, to the random mani- 
fold regime and eventually to the asymptotic Bragg glass 
regime. This nonuniversality could be in principle tested 
by neutron scattering changing the external field. 

We acknowledge discussions with H.E. Brandt, D. 
Feldman and S. Scheidl. This work was supported by 
DFG through SFB 341. 
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